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Determined by Capacitance of Strongly 
Disordered Graphene 

Wei Li''^*, Xiaolong Chen'*, Lin Wang', Yuheng He', Zefei Wu', Yuan Cai^, Mingwei Zhang', 
Yang Wang', Yu Han', Rolf W. Lortz'-^, Zhao-Qing Zhang' ^, Ping Sheng'-^ & Ning Wang'-^ 

^ Department of Physics, ■^the William Mong Institute of Nano Science and Technology, the Hong Kong University of Science and 
Technology, Clear Water Bay, Kowloon, Hong Kong, China. 

We demonstrate that fluctuations of the local density of states (LDOS) in strongly disordered graphene play 
an important role in determining the quantum capacitance of the top-gate graphene devices. Depending on 
the strength of the disorder induced by metal- cluster decoration, the measured quantum capacitance of 
disordered graphene can dramatically decrease in comparison with pristine graphene. This is opposite to the 
common belief that quantum capacitance should increase with disorder. To explain this counterintuitive 
behavior, we present a two-parameter model which incorporates both the non-universal power law behavior 
for the ADOS and a lognormal distribution of LDOS. We find excellent quantitative agreements between the 
model and measured quantum capacitance for three disordered samples in a wide range of Fermi energies. 
Thus, by measuring the quantum capacitance, we can simultaneously determine the ADOS and its 
fluctuations. It is the LDOS fluctuations that cause the dramatic reduction of the quantum capacitance. 



For disordered graphene, theoretical predictions suggest that as the disorder strength increases, the average 
density of states (ADOS) increases accordingly, in comparison with pristine graphene^'^. The changes of the 
ADOS near the charge neutrality point (NP) still remain ambiguous and under debate^"^. One of the 
theoretical predictions on disordered graphene suggested that the ADOS followed a power law, i.e. 
{p{Ep)) ^ \Ep\^ , 0 < a < 1, where Ep is the Fermi energy^'^ and the exponent a decreases with increasing disorder 
. Thus, the ADOS is expected to increase when disorder becomes stronger in graphene. 

In two-dimensional (2D) materials, quantum capacitance^^'^^ can directly reflect the ADOS because it is 
proportional to the ADOS at Ep due to the large contribution of electron compressibility^^, i.e. Cq = e^' (p) , where e 
is the electron charge and (p) is the ADOS at Ep. The measurement of Cq is more immune to the scatterings arising 
from disorder compared to the DC transport measurement. The theoretical expressions of the quantum capa- 
citance and the ADOS in single layer pristine graphene^ ^'^^ have been verified experimentally by capacitance 
measurements for devices with various top-gate geometries 

The measurement of quantum capacitance is expected to help probe the non-conducting states, particularly in 
strong localization systems such as the chemically decorated graphene structures, which contain strong disorder 
as evidenced by the DC transport measurement^^"^^ So far, no experimental work seems to have been reported on 
the quantum capacitance of disorder graphene. Since the ADOS has been predicted to increase in disordered 
graphene, it is natural to expect that quantum capacitance measurements should reflect the increased ADOS 
directly. 

On the other hand, the local density of states (LDOS) in disordered graphene nanostructures has been 
investigated theoretically and experimentally (mostly by scanning tunneling microscopy), and it is believed to 
fluctuate due to the presence of the disorder^^"^^. For substrate-supported pristine graphene, experimental works 
demonstrated that the spatial fluctuations of LDOS in the vicinity of the NP follow Gaussian distributions^°'^\ 
Gaussian distributions have been predicted for weakly disordered systems with extended states^^. However, for 
strongly disordered systems, the fluctuations of LDOS have been predicted to be lognormaP^. Recent numerical 
simulations showed the lognormal distribution in disordered graphene^^. The lognormal distribution has also 
been identified in some experimental studies of 2D materials^^. 

To introduce strong disorder into single layer graphene, we adopted our previously developed method of Ag- 
decoration of graphene surfaces Through transport measurements of Ag- decorated graphene, we found strong 
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Anderson localization induced by the disorder of Ag clusters. The 
resistance of Ag- decorated graphene largely increased right after the 
decoration, and the variable- range-hopping (VRH) model inter- 
preted the experimental data of graphene conductivity fairly well 
even at the temperature as high as 300 K^^. Obviously, strong local- 
ization had been induced in such disordered graphene even at room 
temperature. It should be pointed out that the Ag decoration method 
was proved by the comparison experiments to be non-invasive and 
the lattice of graphene was preserved. Therefore, it makes strongly 
disordered Ag- decorated graphene an interesting candidate for the 
study of both the ADOS and its local fluctuations through capacit- 
ance measurements. 

In this letter, contrary to our expectation, we found that the mea- 
sured quantum capacitance Cq decreases dramatically as disorder 
strength (controlled by the plasma-assisted Ag-cluster decoration^^) 
increases. To understand the reduction in a quantitative way, here we 
proposed a two-parameter model which incorporates both the 
ADOS and its local fluctuations. Based on this model, we determined 
simultaneously the ADOS changes of the Ag- decorated graphene 
and the local fluctuations. It is found that ADOS indeed increase 
with disorder and follow the power law {p{Ep)) ^ \Ep\^.The dramatic 
reduction of the measured Q is caused by large fluctuations of LDOS, 
which follows a lognormal distribution and is independent of the 
Fermi energy. 

Results 

The fabrication of top -gate graphene devices and the measurement 
of capacitance are described in the method section. Fig. 1(a) schem- 
atically shows the geometry of the Ag- decorated graphene device 
with a top-gate. Fig. 1(b) illustrates the three-terminal configuration 
of the capacitance measurement C^( V^c/z) ^nd oxide layer capacitance 
Cox- In the serial connection, the quantum capacitance can be 
obtained from the top -gate capacitance Qg through the relation 
Cq(Vch) = (Ctg(Vtg)~^ — C~^)~^ . Here Vch is the potential across 



the graphene sheet and can be obtained by integrating the charge 
conservation^ relation Ctg(Vtg)'dVtg = Cox'd(Vtg — Vch), i.e., 

Vch = Vtg— —^dVtg. For pristine graphene devices, the 
Jo Cox 2e^kT 
quantum capacitance has been given by CaQ_theory = y 

/ / eVchW ^^^^^^ 

In 2 1+ cosh- — where k denotes the Boltzmann con- 

V V )) 

stant and Vp the Fermi velocity. We fitted the measured quantum 
capacitance of pristine graphene C^o_exp to the theory C^Q theory to 
obtain Cox = 1.14^F/cm^. The result is shown in Fig. 1(c). The con- 
sistency in the vicinity of the NP has been improved compared to 
previous works^^'^^'^°. One of the possible reasons for this improve- 
ment is the use of Cu mask before Yttrium deposition. By this way, we 
effectively reduced contamination from e-beam lithography, which 
may result in electron-hole puddles^°'^\ 

Fig. 1(d) shows the measured quantum capacitances Cq{Vch) of 
three typical Ag-decorated graphene devices. Samples Agl, Ag2 and 
Ag3 were decorated by Ag clusters via sputtering (5 W) for 1 s, 5 s 
and 10 s, respectively. Compared to the results recorded from the 
pristine graphene Qo_exp in Fig. 1 (d), the quantum capacitance 
obtained from Ag-decorated graphene devices is much smaller. 
The Cq{Vch) curves of Ag-decorated graphene devices are all asym- 
metrical between the electron and hole bands. This asymmetry was 
also observed in transport measurements^^. As the Ag decoration is 
increased, the quantum capacitance of the devices decreases dramat- 
ically, which seems to be counterintuitive. Direct interpretation of 
these results via Cq = e^-{p) leads to disagreement with the theor- 
etical predictions of the ADOS in disordered graphene. 

In order to quantitatively understand the dramatic reduction of 
quantum capacitance in disordered graphene samples, we propose 
here a two -parameter model. In this model we incorporate both the 
ADOS of the disordered graphene and its LDOS fluctuations. To 
simulate the quantum capacitance, we first divide the total area A 
of graphene evenly into N small regions A/, each has a size A/N as 




Vch (V) Vch (V) 

Figure 1 | (a) The geometry of a top-gate Ag-decorated graphene device. The small white islands represent the Ag clusters, (b) The equivalent circuit of 
the three-terminal capacitance measurement. The p-Si substrate is grounded to avoid its parasitic capacitance, (c) The quantum capacitance of one 
pristine graphene device with Qx = 1 • 14/xF/ cm^. The inset shows the optical image of the device. The dashed line indicates the outline of the graphene 
flake and the scale bar is 5 [im. (d) The measured quantum capacitance of three Ag-decorated graphene devices sputtered for 1 s, 5 s, and 10 s, 
respectively. 
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shown in Fig. 2(a). For simplicity, we set A = 1. In each region A/, we 
define a normaHzed random variable di = with {Sj) = 1, where 

\Pi/ 

Pi is the LDOS in A/ and (p^ ) is the ADOS of the system. For a given 
Ep, the local quantum capacitance in the region A/ can be written as 

.2 



1 

Cai = e^ O:' — = —•{p;)-di 
^ N N 



(1) 



As indicated by the equivalent circuit shown in Fig. 2(b), the total 
capacitance of the graphene device Q is the summation of all parallel 
capacitances of AT serial capacitances of Q and where Q = CoJN 
is the oxide layer capacitance in each small region. Thus, the 
quantum capacitance can be simulated using the following equation: 



Q_mod ={Ct ^ — (NCo) 



N ' 

Ei/(c,,-+Q) 

i=l 



(2) 



which depends only on the ADOS (pi) and the normalized distri- 
bution of Sf. In order to simulate the quantum capacitance, we need 

to choose a correct distribution for Sj = . In strongly disordered 

\Pi/ 

graphene, the distribution of LDOS has been studied and found to be 
lognormaP^. The lognormal distribution has also been evidenced 
experimentally in other 2D disordered materials^^. Thus, to simulate 
Eq. (2), we propose a lognormal distribution for 3, i.e.. 



m- 



• exp 



2e 



(3) 



or, ^-lniV(-c^V2,^^), [0, + oo) with = Because the 
variance of 3 is exp((^^) — 1, denotes the degree of the fluctuations. 
Theoretically, the ADOS in disordered graphene have been predicted 
to increase with disorder strength. One of the proposed analytical 
expressions for ADOS is a non-universal power law referred to quasi- 
particles in d-wave superconductor, which was supported by numer- 
ical exact diagonaUzations (ED) for different types of disorders^'^. In 
our simulation, we assume that the ADOS {pi) follow the non-uni- 
versal power law. Thus, the ADOS at the Fermi energy {Ep = eVch) is 
given by 




Icjc/- Icjc, 

T T T T 

o 




0.2 0.3 0.4 

V . (V) 



0.5 



Figure 2 | (a) The top-gate Ag-decorated graphene device is divided into N 
small areas in the FIDOS model. Q and C^i are the oxide layer capacitance 
and the quantum capacitance in Aj, respectively, (b) Equivalent circuit of 
the top -gate graphene device in the FLDOS model. The total capacitance is 
the summation of parallel capacitances of all serial capacitances of Q and 
Cqi. (c) The experimental quantum capacitance (scatters) and the 
simulation results by the FLDOS model (solid lines). 



Table 1 


Fitting parameters and DOS at Ep = 


0, 0.25 eV 


Sample 


a 


^ p{0 eV) eV ' m-2 


pyJ.AD By] Gv m 


Pristine ° 


1.00 


0.00 0.00 b 


2.77X 10^^ 


Ag 1 


0.74 


1.66 3.36x10^^ 


7.25 X 10^^ 


Ag2 


0.69 


1.95 9.73 x10^^ 


9.22 X 10^^ 


Ag3 


0.49 


2.75 7.43 x10^^ 


2.41 X 10^^ 


°Only for comparison. '° Vanishes at the NP. 



321/ f 



2e 



1 



(4) 



is the residual ADOS at 



where Pn(o^) = r exp . 

NP, k is the cut-off energy and X = 3t,ais the lattice constant and t is 
the nearest-neighbor hopping energy^'^ {t = 2.8 eV). The exponent 
a G (0,1) in Eq. (4) represents the disorder strength, and its value 
decreases with increasing disorder strength. 

The simulated quantum capacitance Q mod of Eq- (2) now 
depends only on two parameters: the disorder strength a and the 
degree of LDOS fluctuations ^. By using Eq. (2), We performed 
Monte Carlo simulations and found that Q mod always converges 
to a well-defined value if N is sufficiently large for a chosen set of 
(a,0. In Fig. 2(c), we plot the simulated results for all the three 
disordered samples shown in Fig. 1(d). It shows good agreements 
between the simulated and measured quantum capacitance in a wide 
range of Fermi energies. To avoid ambiguity in the vicinity of the 
js^p3o,3i^ the simulation started from Vch = 0.15 V. Good agreements 
are also obtained for the hole branch, (see Supplementary 
Information). However, for Sample Ag3, some deviations between 
the measured and simulated results are noticed. This indicates that 
for a strongly disordered sample, the ADOS might deviate from the 
power law behavior of Eq. (4). 

The fitting parameters for the three samples are listed in Table 1. It 
is clearly seen that a decreases with increasing disorder. The corres- 
ponding ADOS are plotted in Fig. 3(a), which shows explicitly that 
ADOS increases with disorder strength. This is consistent with the 
previous theoretical predictions of ADOS in disordered graphene^'^ 
and our transport measurements^^. Table 1 also indicates that the 
residual DOS po at the NP cannot be ignored in such highly disor- 
dered graphene samples. Thus, the dramatic reduction in measured 
quantum capacitance found in Fig. 1(d) is caused by the large 
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Figure 3 | (a) The simulated ADOS of three samples which follows Eq.(4) 
with a = 0.74,0.69,0.49 for Agl, Ag2 and Ag3, respectively and the ADOS 
of pristine graphene for reference, (b) shows the reduction factor P of 
experimental (scatters) and simulation results by the FLDOS model (solid 
lines). 
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5 = P/<P,> 



Figure 4 | (a) shows the distributions of normalized LDOS at Ep = 0.05 t,0.15 t,0.20 t{t= 2.8 eV), when W = 
distribution, (b) also indicates the single lognormal distribution of LDOS at different £p when W = 1.50 t. 



0.50 ty which all follow one lognormal 



fluctuations of LDOS. Indeed, Table 1 shows that the value of ^ 

increases with disorder strength. 

To quantif)^ the reduction of quantum capacitance, here we define 

C ^ 
a dimensionless reduction factor = — - , where Cq^ = Cqi = 

1 / 
e^-{Pi)- — •^Y2^di = e^-{pi) is the quantum capacitance if the DOS 

is assumed to be uniform. Actually, is also the real quantum 
capacitance of the disordered graphene. As we will show later that 
can only be achieved when the oxide layer capacitance Cox is 
much greater than C^*. In Fig. 3(b), we plot the reduction factor for 
both the measured quantum capacitance, i.e., 



exp 



p/Q^,and 

the simulated quantum capacitance, i.e., i^mod ^ Q-mod/Q*- Here, 
Cox = l'l4:iJ.F /cm^ for all three samples. It is clearly seen that more 
reduction is found for more disordered samples. 

From the above simulation results, we would like to point out 
another interesting property of the LDOS fluctuations in disordered 
graphene. The fact that, for each disordered sample in Fig. 2(c) or 
3(b), the quantum capacitance of disordered graphene can be excel- 
lently described by only two parameters in a wide range of Fermi 
energies indicates that the distribution of the normalized LDOS 3 
given in Eq. (3) is independent of Ep for a given disorder sample. To 
confirm this, we have performed the following simulations based on 
the Anderson tight-binding model for disordered systems^^. We cal- 
culated the LDOS Pi in graphene with a given site disorder W by 
using the Haydock-Heine-Kelly recursion method, (see Supple- 
mentary Information for details). For two different values of W, 
the results of the distributions for three different values of Ep are 
shown in Figs. 4(a) and 4(b), where the hopping energy t is taken as 
the units of energy. These results show that for each given value of 
the distributions of 3 at three different Fermi energies can all be 
described nicely by the same lognormal distribution (solid curve). 
It is also seen that the variance of the distribution is larger for larger 
W. Thus, these simulation results support our two -parameter model 
that the distribution of 3 depends only on the disorder strength and is 
independent of the Fermi energy Ep. 

Finally we would like to study the role played by the Yttrium oxide 
layer in the observed reductions of quantum capacitance of disor- 
dered graphene. The presence of a finite capacitance of the Yttrium 
oxide layer Cox allows local fluctuations of channel potential across 
the graphene sheet in each small region. It is such fluctuations that 
lead to the reduction of the observed quantum capacitance of the 
graphene. For this purpose, we rewrite Eq. (2) in the following form: 

(<5,-/(.5,. + r)> 



where 



In this form, Eq. (5) also exhibits an interesting 



scaling behavior, i.e., the reduction factor depends only on the ratio 
of Cox to Cq*, not on Cox and Cq* individually. Here we focus on one 
particular disorder strength and study the reduction factor as a func- 
tion of r by using different thicknesses of the Yttrium oxide layer. We 
also test the r- scaling law by using data obtained at different Fermi 
energies, where Cq* varies. 

We measured the quantum capacitance of the same Ag- decorated 
graphene from devices with different Cox t>y tuning the thickness of 
the Yttrium oxide layer. The decorations of all these samples were 
carried out simultaneously under the same conditions as Sample Agl 
in order to guarantee the same strength of disorder. The results of 
reduction factor for different Yttrium oxide layer thicknesses, which 
correspond to = 1.13, 0.89, 0.76, 0.70^F/cm^ and different 
Fermi energies {Vch = 0.15, 0.25, 0.35, 0.45 V) are shown in Fig. 5. 
All these data overlap nicely with the simulated solid curve obtained 
by using a = 0.74, ^ = 1.66 shown in Table 1 for Sample Agl. Fig. 5 
verifies the r-scaling behavior of our two-parameter model with a 
lognormal distribution of LDOS. We also notice that the experi- 
mental data slightly deviates from the model when Cox becomes 



0.24 




0.06 



0.00 



0.04 0.06 0.08 0.10 0.12 0.14 0.16 

r = C /C . 

ox q* 

Figure 5 | The inset shows the reduction factor P (calculated by the FLDOS 
model) reaches 1 when Cg^^s very large at different ^. The reduction factor 
P (scatters) of graphene decorated by Ag for 1 s at 5 W with 
Cox = 1.13, 0.89, 0.76, OJOjuF/cm^ agrees with the r-scaling of the FLDOS 
model (solid lines). 
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smaller. In this case, the oxide layer becomes thicker and the induced 
polarization effect might enhance the fluctuations of LDOS, leading 
to a smaller In the inset of Fig. 5, we show the function obtained 
from different values of c^. It is found that decreases with increasing 
fluctuations of LDOS as expected. It also shows that p increases with r 
and saturates to 1 at very large r (unachievable practically), indicating 
that the real quantum capacitance Cq^=e^-{p) of disordered gra- 
phene could only be achieved when Cox ^ Cq*. 

Discussion 

Recently, experimental works demonstrated that the spatial LDOS 
fluctuations in the vicinity of the NP follow a Gaussian distri- 
bution^°'^\ Compared to the lognormal distribution which can 
describe LDOS in strongly disordered systems, the Gaussian dis- 
tribution is only capable of describing systems with weak disorder. 
We have also performed simulations of Eq. (5) based on the 
Gaussian distribution. However, we do not found any significant 
reduction in quantum capacitance, (see Supplementary Informa- 
tion). In strongly disordered systems, the envelope of the localized 
wave functions decays exponentially from the center of the local- 
ization sites^^'^^, creating strong spatial fluctuations of LDOS. This 
results in a dominance of LDOS at lower values than its mean and 
gives rise to a positive skewness. Although the variance of both 
Gaussian and lognormal distributions can represent the degree of 
fluctuations, the positive skewness of the lognormal distribution 
could be the underlying reason why it can correctly describe the 
LDOS fluctuations in disordered graphene, where significant 
LDOS fluctuations occur. 

In conclusion, we have established a simple two -parameter ( a, 0 
model for the quantum capacitance of disordered graphene. This 
model incorporates both the power-law dependence of the ADOS 
and the lognormal distribution of LDOS. By measuring the quantum 
capacitance of disordered graphene we simultaneously determine 
both the ADOS and its local fluctuations of the system. The model 
gives quantitative agreements with the measured quantum capacit- 
ance of three different strongly disordered Ag- decorated graphene 
samples. For Ag-decorated graphene samples, a given Ag-decoration 
corresponds to a configuration of disorder induced in the graphene 
sheet. The presence of disorder increases the ADOS and also induces 
large spatial fluctuations in LDOS compared to pristine graphene. 
The reduction of measured quantum capacitance is caused by large 
fluctuations of LDOS mediated by the presence of the insulating 
oxide layer. The validity of the model is also supported by our simu- 
lation results of LDOS using the Anderson model that the distri- 
bution of normalized LDOS is lognormal and independent of 
Fermi energy for a fixed disorder. 

Methods 

Device fabrication. The pristine graphene devices were prepared on p-Si 
substrates coated with 300 nm Si02 insulating layers, as reported previously^**. The 
top-gate electrodes were then added using a standard e-beam lithography 
technique^^'^**'^''. The insulating layer between the top-gate electrode and graphene 
was fabricated by depositing a thin layer of Yttrium on the graphene devices, 
followed by oxidation at 180°C for 10 min in air. We used a mask (thin Cu foils) 
in order to minimize the influence of direct electron beam irradiation^^ when 
patterning the insulating layer. The silver clusters were first deposited on the 
pristine graphene by the method reported previously^**. Then, the graphene 
devices were immediately transferred into a high vacuum chamber for Yttrium 
deposition by electron-beam evaporation. 

Capacitance measurement. The capacitance was measured by a Keithley CV 
analyzer (model 590) with a sensitivity of ~0.1/F. As illustrated in Fig. 1(b), the 
capacitance measurement was performed in a three-terminal configuration and 
the p-Si substrate was grounded to avoid its parasitic capacitance. The residual 
capacitance of the circuit in the same geometry of the devices without graphene is 
smaller than 0.3/F. The area of graphene for capacitance measurement was 
normally larger than 30 /im^ as shown in the inset in Fig. 1(c). Pristine graphene 
devices were used to calibrate the capacitance of the Yttrium oxide layer (QJ of 
the Ag-decorated ones in the same batch. We extracted the quantum capacitance 



of a large number of devices for the same Ag decoration in different batches, and 
the results were similar. 
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